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Abstract 

A brief discussion of massive tadpole diagrams and their phenomenological consequences is presented. This in¬ 
cludes predictions of the p parameter and, as a consequence, the mass of the W boson, implications on the charm and 
bottom quark masses from the moments, i.e. the derivatives of the current correlators, and the Higgs boson decay 
rate. A fairly consistent picture emerges, with m c (3 GeV) = 0.986 ± 0.013 GeV and = 3.610 ± 0.016 GeV. 

Furthermore, fairly stringent predictions for the Higgs decay rate into photons and gluons are obtained, which will be 
interesting in increasingly precise experiments. 
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1. Introduction 

Using massive tadpole diagrams significant progress 
has been made in the improved prediction of various 
physical quantities. They have, e.g., contributed to an 
amazingly precise relation between quark masses (in 
particular the top-quark mass M t ), the mass of the Higgs 
boson, Mh, and the masses of the gauge bosons, M w 
and M z . Three- and even four-loop corrections have be¬ 
come accessible during the past years. Direct and indir¬ 
ect measurements are well consistent, at least within the 
current world average M t = 173.21+0.51+0.71 GeV HI 
and Mh - 125.7 + 0.4 GeV (T]|- Current-current cor¬ 
relators, evaluated in three- and partially four-loop ap¬ 
proximation are thus an indispensable tool for tests of 
the Standard Model (SM), as we will discuss in more 
detail in Section [2] 

The evaluation of three- and even four-loop tad¬ 
pole diagrams is directly related to the evaluation of 
moments of the charm- and bottom-quark correlators. 
These may in turn lead to a precise determination of 


* Corresponding author 

Email address: matthias. steinhauser@kit. edu (Matthias 
Steinhauser) 


charm- and bottom-quark masses. Since all these quant¬ 
ities, in turn, are directly accessible, both to a per¬ 
turbative and a non-perturbative treatment, a remark¬ 
ably consistent picture emerges. In fact the analysis for 
the charm- as well as the bottom-quark mass leads to 
m c (3 GeV) = 0.986 + 0.013 GeV and m*(10 GeV) = 
3.610 + 0.016 GeV El El, a result quite comparable to 
other methods, in particular to non-perturbative studies 
as will be discussed in Section [3] 

Finally, in Section[4| we list a collection of topics and 
results which are connected to the main theme of this 
article. This includes the decoupling of heavy quarks 
at four loops, the Higgs-gluon coupling up to five-loop 
order, and the Higgs-decay rate into two photons at three 
and four loops, including non-singlet and singlet terms. 

2. Weak corrections 

Let us start with the p parameter, the quantity intro¬ 
duced in the early times of electroweak interactions of 
the SM. In its more modern version it gives the relation 
between gauge boson masses, the weak- and the electro¬ 
magnetic coupling Gp and a , a relation, which exists at 
tree level already. In higher orders this relation depends 
on all remaining parameters of the SM. The radiative 
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corrections are dominated by the quadratic dependence 
on the mass of the top quark M t , the logarithmic de¬ 
pendence on the Higgs boson, M H , and, to a lesser ex¬ 
tent, also on the masses of the remaining quarks m q and 
leptons mf. 

Mw = f(G F , M z , or, M u M H ; m q , mi). (1) 


three-loop regime was taken in the limit M H = Ofl.In 
fact, this turns out to be a poor approximation, leading 
to tiny corrections for the terms of order Xf and a s Xf 
with 


G f M] 
8 V27T 2 


( 5 ) 


A slightly different version of the same equation 


M\ 





na 

V 2G F 


(1 + Ar) 


( 2 ) 


makes the presence of the electroweak corrections even 
more transparent. This equation can be rewritten and 
simplified even further by separating A r into a piece 
which is dominated by weak effects and another one 
which is dominated by electromagnetic effects, mainly 
due to the running of the electromagnetic coupling A a 
(see Refs. 0 and 0 for three- and four-loop correc¬ 
tions, respectively). Furthermore, it is convenient to 
separate the leading M] dependence which leads to 


COS 2 0\y 

A r = —— r -A p + A a + A r remlillillg . 

sin 6 W 


( 3 ) 


The first three-loop result with M H different from zero 
requires the full set of three-loop tadpoles Go). At order 
a s Xjf(M t /M H ) this corresponds to QCD corrections to 
the two-loop diagrams of order Xjf(M t /M H ) (see Fig. [l] 
for a sample Feynman diagram). At order X^ f(M t /M H ) 
diagrams with one quark line contribute, as well as those 
involving two disconnected quark lines. 



Figure 1: Sample diagram for the a s Xff(M t /MH ) contribution. Solid 
lines denote top quarks, dashed lines Higgs or Goldstone bosons and 
curly lines gluons. 


Here A p = 3G F M?/(8 V27r 2 ) + ... incorporates the dom¬ 
inant weak terms evaluated in leading order by Velt- 
man [61 nearly 40 years ago. 

It is the aim of ongoing and of future theoretical stud¬ 
ies to compete with the precision anticipated for the next 
round of experiments. The present precision and the 
precision anticipated for the future (as far as SM W and 
SM t are concerned) are given by 


5M W [MeV] 

6M t [GeV] 


33 

5 

status 2003 (LEP, TEY) 

15 

0.76 

now (TEVATRON, LHC) 

8—>5 

0.6 

aim (LHC); theory 

3—>1.2 

0.1-0.2 

ILC, TLEP 


As it turns out, the relative shifts of M w and M t are just 
of the right order of magnitude to explore the sensitivity 
towards radiative corrections. This is seen most clearly 
by considering the shift in M t that is compensated by a 
shift in M w 

6M W * 6 • 10 - 3 8M, , (4) 

keeping a(M z ), M z and M H fixed. 

Let us now recall the development of the theory pre¬ 
dictions for A p during the past one or two decades. 

Early results related to the two-loop approximation 
can be found in Refs. dill- These papers are based on 
the approximation Mf » M 2 W . The first step into the 


At the same time the translation from the MS mass 
m t (M t ) to the pole mass M t has to be performed at two 
loops. This corresponds to the evaluation of two-loop 
on-shell diagrams. For the case of the Xff(M t /M H ) cor¬ 
rections the counterterms are of order Xj and are depic¬ 
ted in Fig. [2] 



__ * 

H, X 


Figure 2: Two-loop on-shell diagrams which contribute to the transla¬ 
tion of the MS mass to the pole mass. 

In contrast to the pure QCD problem (see Section |3J 
with only one mass scale being present (and which can 
therefore be expressed in closed analytic form) here 
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one typically encounters two different scales, M H and 
M t . A closed analytical solution is no longer at hand. 
There are, however, various cases which allow for ex¬ 
pansions, one of which is perfectly valid for M H and M t 
in the region of interest. Analytic results are available 
in the cases M H = 0 and M H = M u where only one 
scale is present. Expansions, which in principle allow 
for arbitrary high precision are then accessible in two 
cases: for the case of large M H with the approximation 
in (M 2 / Mjj) n modulo logarithms which is valid down 
to M h « 2 M t and the case of M H around M t which 
is valid from fairly small M H , say M H « 0.1 M u up to 
Mh ~ 2M r . The results for the expansion in (Mh / M t ) n 
and in (M H - M t ) 2 / M 2 are shown in Fig. [ 3 ] Note that for 
Mh = 0/126 GeV one obtains for the prefactor of a s X 2 , 
a part of the three-loop term of Ap, the values 2.9/120. 


•••/ 


Ap' 


up to 1. ord (m ~M n ) 
up to 2. ord (m ~M n ) 
up to 3. ord (m t ~M H ) 
up to 4. ord (m ~M n ) 
up to 5. ord (m t =M H ) 
0. ord (m t <M H ) 
up to 1. ord (m t <M H ) 
up to 3. ord (m ( <M H ) 
up to 5. ord (m t <M H ) 


(a s /7t)X"Ap 


Figure 3: Contributions of order a s Xj to Ap in the on-shell definition 
of the top-quark mass. The black squares indicate the points where 
the exact result is known. Using the latest numerical values for M t 
and Mh one obtains Mn/M t ~ 0.73. 

The results for the shift in M w and the effective weak 
mixing angle are shown in Fig. [4] for the four contribu¬ 
tions which are most relevant. The terms proportional to 
X t and X t a s must be taken into account in any sensible 
analysis and amount to a shift in dp of order 0.00865. 

The two-loop piece proportional to X 2 is of the same 
order as the three loop piece, proportional to a 2 X t . The 
purely weak term proportional Xf is negligible now and 
in the foreseeable future, the term proportional to a s X 2 
is just below the present sensitivity. 

Four-loop QCD contributions 

Two- and three-loop QCD corrections to Ap have 
been computed in Refs. EH about 20 years ago. 
As stated above, since several years it is now possible 
to push the predictions for the p parameter to the four- 
loop level. This requires, on the one hand, the relation 
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Figure 4: The shift in Mw and sin 2 0 e ff as a function of Mn/M t 
induced by the corrections of order Xf, a 2 s X t , a s X t , and Xj. 


between pole- and MS mass in three-loop approxima¬ 
tion ED HE) and, on the other hand, the evaluation of 
about fifty four-loop tadpole diagrams. This has been 
achieved by a combination of analytically methods, dif¬ 
ference equations and semi-numerical integration [ 1 19;- 
21]. In a first step this has lead to the four-loop result 
for the p parameter in the MS scheme 

( 4 ‘oops) = 3 Gf " 1 > ( M 3 (-3 2866 + 1.6067) . ( 6 ) 

8 yfljr 1 V > ' ---' 

- 1.6799 

The first term has been evaluated in Ref. l22l . the 
second one, which (in the MS scheme) leads to a re¬ 
duction by a factor of about 1/2, in Refs. 119, 20]. For 
the numerical evaluation the translation from the MS to 
the pole mass is more convenient and one finds 

6 p‘ 4l0ops) = 3 9zMl. (2A 3 (-93.1501) (7) 

8 V27T 2 V a -) 

which corresponds to a shift in M w by about 2 MeV, 
similar to the three loop term of order a s X 2 . 


3. Charm- and bottom-quark masses 

The precise determination of the charm- and bottom- 
quark masses from relativistic four-loop moments can 
be considered as one of the truly remarkable successes 
of quantum field theory with remarkable agreement 
between perturbative EH and lattice methods [23 - 25 ]. 
Fet us first give a more detailed motivation, then present 
the theoretical framework and finally compare the res¬ 
ults based on perturbative and lattice methods. 

Precise bottom quark masses enter many B physics 
quantities. During the past years significant progress 
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has been made on the one hand in the analysis of B- 
meson decays and, on the other hand, the T spectro¬ 
scopy. In particular the latter has led to a fairly consist¬ 
ent result of m b (m b ) = 4.193^ ^ GeV (26) (see also 
Ref. (27) where m b (m b ) = 4.169 ± 0.009 GeV has been 
obtained) in excellent agreement with the result based 
on sum rules discussed in detail below. 

Let us further motivate the need for precise quark 
masses for the case of the bottom-quark mass, which 
enters a number of physical observables. Most promin¬ 
ently, we want to mention the decay of a Higgs boson 
into bottom quarks which, using the scalar correlator to 
five-loop precision (28), can be written in the form 


T(H -> bb) 
R (S \M h ) 


GrMj, 

4 yfln 
1 + 5.667 


+ 41.758 


m 2 b (M H )R (s \M H ) ( 8 ) 

( 7 )—(?) 2 


825.7 

1 +0.19551 +0.03469 
+ 0.00171 -0.00117. 


(9) 


The theory uncertainty, which is generously taken from 
a variation of the scale parameter between M H /3 and 
3 Mh, is reduced from 5%o for the four-loop to 1.5%o 
for the five-loop result. Thus, the main uncertainty is 
induced from the uncertainty of the bottom quark mass 
which (at energy scale 10 GeV) is given by 0 

/ ot c 0.1189 \ 

m^lOGeV) = 3610- —-12+11 MeV. 

v ' \ 0.002 / 

The running from 10 GeV to M H depends on the an¬ 
omalous mass dimension y m , the /3 function and on a s . 
With the present knowledge (i.e. four-loop anomalous 
dimensions as implemented in RunDec 129,, !30:|) one 
finds 


m b (M H ) = 2759 + 8 | m , + 27| ffj MeV. (10) 


It is interesting to investigate the effect of the five- 
loop anomalous dimensions. Taking into account the 
y m to five-loop accuracy [31| 32 ] together with ^ 4 , the 
five-loop contribution to the function of QCD, which 
is still unknown, one obtains the following uncertainties 


4m 2 h (M H ) 

m 2 b (M H ) 


-1.3 X 10- 4 OS 4 /A) = 0) 

-4.3 x 10" 4 OS 4 /ySo = 100) 

-7.3 x 10 -4 (/? 4 //?o = 200), (11) 


which lead to an uncertainty of a few MeV in m b (M H ) 
which is small compared to the current error shown in 
Eq. ^ . 


Another motivation which also points to a precision 
around 10 MeV is based on the picture of Yukawa uni¬ 
fication. In this approach A r ~ A b ~ A t at the GUT 
scale. In effect this implies ~ Sm t /m t . Assum¬ 

ing a precision of the top-quark mass Sm t ^0.5 GeV 
then leads to a precision of Sm b ~ 10 MeV, consistent 
with our finding below. 


SVZ sum rules, moments and tadpoles 

The main idea originally advocated in Ref. (33) is 
based on the observation (cf. Fig. [5]), that the cross 
section for hidden (//'F plus V F(2S')) plus open (DD 
plus resonances) charm production is well described by 
perturbation theory, if the analysis is restricted to suffi¬ 
ciently low moments. 
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Figure 5: Sketch of the R ratio in the charm-threshold region. 


Let us first recall some definitions. The two-point 
correlation function 

(~<l 2 gnv + <Mv)n(<? 2 ) = i J' dxe !qx (0\T j M (x)j v (0)\0) 
is related to the electromagnetic current j^ as follows 
R(s) = \2n\m [ll(q 2 = s + ie)] . (12) 

In fact, we are only interested in lowest moments of n, 
corresponding to the first few terms of the Taylor expan¬ 
sion of n (q 2 ): 


n< <r ) = e 5+2 c "‘‘- <>3) 

n> 0 

where Q q corresponds to the charge of the considered 
quark. Here z = q 2 / (Am 2 q ) and m q = m q (p = m q ) is the 
MS mass at the scale p = m q . Let us, for definiteness, 
restrict the following discussion to the charm quark, i.e., 
q — c. 
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For the moments one finds 
C n = Cf +^cl 1) + (^) 2 cf ) 

/0's \ 3 t^(3) 

+ ( 7 )c. +.... 

if the renormalization scale is set to /i = m q or 


(14) 



if one is interested in the generic form with 
l mc = In(m 2 (yu)/yu 2 ). The next-to-next-to-leading order 
calculation had been performed already nearly twenty 
years ago 134 - 361 and is available for all four (vector, 
axial, scalar and pseudoscalar) correlators. The original 
evaluation was up to n = 8 . More recently, this has 
been extended to n = 30 [37, 38]. Now this project has 
been pushed to N 3 LO, and the results will be described 
in the following. In a first step the nj. contribution has 

been computed for Co and C\ [39], then the complete 
result became available. The reduction of the many dif¬ 
ferent diagrams has been performed to 13 master integ¬ 
rals, shown in Fig. [ 6 } using the Laporta algorithm (40). 
Subsequently these 13 remaining integrals are evalu¬ 
ated, using originally a combination of numerical and 
analytical, now purely analytical methods. 



Figure 6: The 13 master integrals. Solid lines denote massive propag¬ 
ators; dashed lines represent massless propagators. 


n 

M r w es 
xl0 (n 1) 

^jthresh 

xl0 (n 1) 

yVfonl 

xl0 (n 1) 

Air 

xl0 (n 1) 

AC 

xl0 (n 1) 

1 

0.1201(25) 

0.0318(15) 

0.0646(11) 

0.2166(31) 

-0.0001(2) 

2 

0.1176(25) 

0.0178(8) 

0.0144(3) 

0.1497(27) 

0.0000(0) 

3 

0.1169(26) 

0.0101(5) 

0.0042(1) 

0.1312(27) 

0.0007(14) 

4 

0.1177(27) 

0.0058(3) 

0.0014(0) 

0.1249(27) 

0.0027(54) 


Table 1: The experimental moments in (GeV) -2n as defined in 
Eq. fl~8) are shown, separated according to the contributions from the 
narrow resonances, the charm threshold region and the continuum re¬ 
gion above - 4.8 GeV. The last column gives the contribution 
from the gluon condensate. 


The reduction of hundreds of integrals to master in¬ 
tegrals has been achieved, originally for Co and C\ RT1 
421, subsequently for C 2 and C 3 , using the program 
Crusher [43] 44]. In the meantime all master integ¬ 
rals are known in closed analytic form to high order in 
s, using results by a number of different authors (see 
Refs. USED and references therein). The results for 
C 4 up to C 10 are known approximately, using an ap¬ 
proximation based on additional information from low 
energies (< q 2 = 0), from threshold ( q 2 = 4m 2 ) and 
from high energy 14914511 . Closely related results are 
also known for axial, scalar and pseudoscalar correlat¬ 
ors [0123E2I- These can be used for the investigation 
of correlators on the lattice f23l and will not be investig¬ 
ated further in this more phenomenological study. The 
heavy quark vector current correlator for q 2 <^c m 2 has 
also been determined in the large limit ED in order 
to study the large-order behaviour. 

The moments are directly related to measurements as 
follows. From theoretical considerations one finds 


At 


127T 2 

n\ 




(16) 


where the quantity C n depends on a s (ji 2 ) and In(m 2 /yu 2 ). 
As default value for yu we use /i = 3 GeV. 

To obtain experimental moments one considers the 
correlator given by 


9 q 2 r R c (s ) 

U c (q 2 ) = -^3 J clv-^——+ subtraction, (17) 


which leads to 


M 


exp 

n 


f 


o «+1 CV 7 


(IB) 


Imposing the constraint At P = At leads to m c at the 
scale fi = 3 GeV, a result that could be easily translated 
to arbitrary values of /i. 
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Figure 7: The normalized cross section R(s ) between 2 GeV and 
10 GeV. The solid line corresponds to the theoretical prediction. The 
uncertainties obtained from the variation of the input parameters and 
of ju are indicated by the dashed curves. The inner and outer error bars 
give the statistical and systematical uncertainty, respectively. The data 
points are from BES |54 55], MD-1 ED and CLEO gT). The vertical 
dashed lines correspond to the location of the J/Y and resonances. 


Let us first discuss the ingredients for charm, and 
then for bottom quarks. The results for the electronic 
widths of J/Y and Y' are taken from the combination 
of BES, CLEO and BABAR experiments, and for the 
continuum R(s) from BES. For the charm case there is 
also a non-perturbative contribution which is, however, 
negligible for the three lowest moments and remains re¬ 
latively small even for the fourth. A careful investiga¬ 
tion of non-perturbative terms combined with the extra¬ 
polation of R u ds as well as R c in the region sufficiently 
far above the respective threshold leads to a remarkably 
consistent result, with errors on m c (3 GeV), as extrac¬ 
ted for n = 1 to 3, below 20 MeV. The result for the 
moments individually split into contributions from the 
narrow resonances, threshold and continuum region is 
shown in Tab. |T| R(s) around the charm threshold re¬ 
gion is shown in Fig [7] 

In particular we observe a remarkable consistency 
between the results for n = 1,2,3 and 4 and a relat¬ 
ively small shift when moving from three- to four-loop 
approximation (cf. Fig. [8]). 


n 

m c (3 GeV) 

exp 

a s 


np 

total 

1 

986 

9 

9 

2 

1 

13 

2 

976 

6 

14 

5 

0 

16 

3 

978 

5 

15 

7 

2 

17 

4 

1004 

3 

9 

31 

7 

33 


Table 2: The second column shows the results for m c {3 GeV) in MeV. 
The errors in the four inner columns are from experiment, a s , vari¬ 
ation of n and the gluon condensate. The last column shows the total 
error. 


In taking the lowest moment as our final result we 
find (3) 


(19) 


When converted from // = 3 GeV to the scale m c , this is 
modified to m c (m c ) = 1279+ 13 MeV, nicely consistent 
with other determinations l23l 24]. The robustness of 
our result is demonstrated in Fig. [8] where the results are 
compared for different orders 0(a l s ), with i-0 , 1, 2 and 
3 and for different moments, with n varying between 
n- 1 and n — 4. 



Figure 8: Dependence of m c (3 GeV) on the number of moments n and 
on O for i = 0,..., 3. 


The result can be compared to those from a large 
number of results, which are based on various differ¬ 
ent observables. Fig [9] shows a compilation of recent 
analyses • 



Figure 9: Comparison of m c {3 GeV) with several other results. 


Similar considerations are applicable for the cor¬ 
responding investigations of the T-resonances and the 
mass of the bottom quark. For convenience of the 
reader we again list in Tab. [3] separately the contribu¬ 
tions from the narrow resonances (Y(1S) - T(45')) (67), 


m c (3 GeV) = 986+ 13 MeV. 
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the threshold region (10.618 GeV-11.2 GeV) (68) and 
the perturbative continuum (E > 11.2 GeV). 


n 

^res,(IS-4S) 

xio (2n+1) 

^thresh 

xiqO+O 

Ajcont 

x 10 (2b+1) 

M exp 

x 10 (2b+1) 

1 

1.394(23) 

0.287(12) 

2.911(18) 

4.592(31) 

2 

1.459(23) 

0.240(10) 

1.173(11) 

2.872(28) 

3 

1.538(24) 

0.200(8) 

0.624(7) 

2.362(26) 

4 

1.630(25) 

0.168(7) 

0.372(5) 

2.170(26) 


Table 3: Moments for the bottom quark system in (GeV) 2n 


For the lowest moment the latter gives the main con¬ 
tribution, starting from the second moment the reson¬ 
ance and the threshold regions are again dominant. In 
particular moments number two and three offer a fair 
compromise between smallness of the error and the 
contribution of the threshold region. Significant pro¬ 
gress has been made between the first measurements 
of CLEO (69) and the more recent one by the BABAR 
collaboration in particular in the continuum region. As 
expected 0 the original CLEO result is too large by a 
factor around 1.3 but reproduces well the qualitative be¬ 
haviour. The more recent BABAR result (68), however, 
is significantly more precise and well in agreement with 
the expectations for the continuum, based on the parton 
cross section. Let us note in this connection that the ori¬ 
ginal BABAR result for R b (s) has to be deconvoluted 
with respect to initial state radiation (ISR), a fact that 
leads to a slight shift of R b (s). 

The results for the bottom quark mass for the lowest 
four moments are given in Tab. [4] (3] IbTTl . For our final 
results we choose n- 2 which leads to 

m b (m b ) = 4.163 ±0.016 GeV, 

m*(10 GeV) = 3.610 ± 0.016 GeV, 

m b (M H ) = 2.759 ± 0.028 GeV. (20) 


n 

m*(10 GeV) 

exp 

<*s 

L* 

total 

m b (m b ) 

1 

3597 

14 

1 

2 

16 

4151 

2 

3610 

10 

12 

3 

16 

4163 

3 

3619 

8 

14 

6 

18 

4172 

4 

3631 

6 

15 

20 

26 

4183 


Table 4: The different columns show the results for m b { 10 GeV) in the 
second column, obtained from the different moments listed in the first 
column. The last column gives the value of m b (jn b ). The three inner 
columns give the uncertainty due to the error in the experimental mo¬ 
ments (exp), the uncertainty due to the error in a s and the uncertainty 
due to the residual scale dependence p. The second to last column 
gives the total uncertainty. All masses and uncertainties are in units of 
MeV. 

The consistency, when comparing the results for the 
lowest four moments is very close to the one from 


2007 0, where only estimates were available for the 
four-loop term of n = 2, 3 and 4. Furthermore, only 
recalibrated results for the continuum corresponding to 
the aforementioned factor 1.3 were available. The res¬ 
ult for mbirrib) can also be compared to those from other 
studies in Fig. [10] OHH ESEl |58l El IM1 170H731 . Al- 
though somewhat towards the low side, the result are 
well compatible with those of earlier investigations. 



m b (m b ) (GeV) 


Figure 10: Comparison of m^(m^) with several other determinations. 

In Fig. |TT] the results for the charm and bottom 
quark masses as obtained from the low-moment sum 
rules (2j|3 | 58, 611 are compared with the numerical val¬ 
ues proposed by the PDG for the years between 2000 
and 2014. It is interesting to note that the extracted 
mass values, which were first based on three, later on 
four-loop perturbative input, remained rather constant 
whereas the PDG numbers seem to converge towards 
these results. 

4. Further applications of massive tadpoles 

Decoupling function at four loops 

In many QCD applications the mass of a heavy quark 
m is much larger than the characteristic momentum 
scale y/s of a considered physical process. As a res¬ 
ult these different mass scales involved in the process 
can lead to potentially large logarithms like lo g(p/ yfs) 
or lo g(p/m) when using an MS-like renormalization 
scheme. In such a situation one can not set the renor¬ 
malization scale p to two different mass scales simul¬ 
taneously, so that a proper choice of p in order to avoid 
large logarithms is not possible anymore. However, by 
“integrating out” the heavy quark field one can construct 
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Figure 11: Comparison of m c (m c ) and as obtained using low- 

moment sum rules Emms] ed (narrow band) and the values from 
the PDG between 2000 and 2014. 


an effective field theory with ni = rif — 1 light quarks 
only, where rif is the number of quark flavours. 

The MS coupling constants and of the 
quark-gluon interaction in the full ^/-flavor QCD and 
the effective ^/-flavor one are different and are related 
by the decoupling function £ g (p, (ji), ni) through the 

matching condition 


for reads 1761 



where 

4’ = i + <(0), 

^ = i + n“(0), 

I? = l + r Gcc(0,0) 


( 22 ) 


(23) 


with the gluon G and ghost c vacuum polarization 
functions U^(q 2 ) and Il® h (q 2 ). The vertex function 
T^ciP, k) is the one-particle irreducible part of the am¬ 
putated Green’s function, where p and k are the outgo¬ 
ing four momenta of the fields c and G, respectively. 
The computation of these functions leads again to the 
evaluation of tadpole diagrams up to four-loop order. 
The four-loop contribution can be expressed in terms 
of the 13 master integrals shown in Fig. [6] The renor¬ 
malized decoupling function is obtained from 


a 


(m) 


( y( n f ) 

A? 


7 
V A? 


/° 
im) g 


Os f) (r) = t? g a ( s f) in). 


(24) 


where (Z g ) 2 is the renormalization constant of the strong 
coupling constant. 

The decoupling function plays an important role in 
testing QCD by running the strong coupling to differ¬ 
ent energy scales. For example, the strong coupling 
a s {m T ) can be measured at the scale of the r-lepton 
mass m T . In the next step one can run this value to the 
scale of the Z-boson mass M z by using the proper run¬ 
ning and decoupling at the heavy charm- and bottom- 
quark thresholds and compare it to the experimentally 
measured result of a s (M z ). This procedure provides 
thus an excellent test of QCD asymptotic freedom. 
The four-loop contribution to the decoupling function 
leads to a reduction of the matching-related uncertain¬ 
ties in the evolution of a s (m T ) to a s (M z ) by a factor of 
two mm. 




(21) Higgs-gluon coupling to N 4 LO 


At leading order the decoupling function is equal to 
one, but receives corrections in higher orders of per¬ 
turbation theory. This matching condition for the MS 
strong coupling constant a s at a heavy quark threshold 
has been computed in Refs. [74, 75] to four-loop order. 
The decoupling function can be determined through the 
computation of polarization functions. The bare relation 


Gluon fusion is the dominant production mechanism 
of the SM Higgs boson H at the Large Hadron Collider 
(LHC), where the leading order process is already at the 
one-loop level and the Higgs boson is produced by the 
fusion of two gluons through a heavy top-quark loop. 
The decoupling function enters in this context as an im¬ 
portant building block since it can be used to derive the 
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effective coupling of the Higgs boson to gluons via the 
following low-energy theorem 


1 2 1 dm° 

Ci enters into an effective Lagrangian 

Xeff = -2 1/4 G l J 2 HCdO[\ 


(25) 


(26) 


in QCD with five flavours, where the top mass depend¬ 
ence is contained in Ci. The symbol Gp is the Fermi 
constant, and [Gj] is the renormalized form of the op¬ 
erator 0\ - G^ V G^ V , where G^ y is the gluon field 
strength tensor. The prime indicates that the object is 
in the effective five-flavour theory and the superscript 0 
denotes a bare quantity. Using the four-loop result for 
f 2 of Refs. (74] 75) allows to determine Ci in four-loop 
approximation, which confirms the result of Ref. m 
in a completely different and independent way. With 
the help of the anomalous dimensions even the five-loop 
contribution to C\ has been predicted [74. 751 up to un¬ 
known five-loop ^/-dependent terms of the p function. 


Decoupling of heavy quarks from the running of the fine 
structure constant 

In complete analogy one can determine from the 
massive photon vacuum polarization function the 

photon decoupling function 

~ i + n“(0)' (27) 

The three-loop results of Ref. l76l have been extended 
to four loops in Ref. [771. Starting from three-loop or¬ 
der there arise also diagrams where the external photon 
couples to massless fermions with the insertion of a 
heavy fermion loop. At four-loop order also singlet 
type diagrams arise for the first time, where the photon 
couples to two different fermion loops. Some example 
diagrams are shown in Fig. [12] 

Higgs boson decay to photons 

The photon decoupling function can be used to de¬ 
termine higher order QCD corrections to the partial de¬ 
cay width of the Higgs boson into two photons (y). The 
amplitude of the partial decay width H —> yy 

T(H -> yy) = ^ \a w (t w ) + Y J A f {r f ) f , (28) 

consists of two parts, a purely bosonic part A w (tw) and 
a purely fermionic part Af(rf) with t w = M^/(4M^) 




Figure 12: The first two diagrams are example singlet diagrams; 
the last four diagrams are examples for the situation where the ex¬ 
ternal photon couples to a massless fermion loop with the insertion of 
a heavy internal fermion loop. The solid lines represent heavy top- 
quarks, the twisted lines denote gluons, and the dashed lines represent 
massless quarks. 


and Tf- M 2 H /(4M 2 ). Higher order QCD corrections 
modify the fermionic part Ay( 17 ) of the amplitude. The 
top quark gives a dominant contribution to the amp¬ 
litude Af (/ = t), since it is the heaviest fermion in the 
SM. In the heavy top-quark mass limit one can again 
describe the Higgs-photon-photon interaction in terms 
of an effective Lagrangian approach 

zjO 

£ eS = —^C° ly F'°^F'°, (29) 

with the vacuum expectation value v° and the field 
strength tensor F'® v . The subscript 0 indicates a bare 
quantity and the prime denotes that the quantity is con¬ 
sidered in the effective theory with ni light active quark 
flavours. The coefficient function depends on the 
photon decoupling function 


C 


o _ 
iy “ 



(30) 


This approach allows one to determine the leading con¬ 
tributions in the heavy top-quark mass limit to the 
Higgs-boson decay into two photons, where the external 
photons couple to the same heavy fermion loop, the so 
called non-singlet contributions. 

At three-loop order in perturbative QCD the non¬ 
singlet contributions to the decay H —> yy have been 
computed in Ref. (781 with several different methods, 
including power corrections of the order [M^/(4M 2 )] 2 . 
The singlet contributions have been added in Ref. (751, 
where also additional power corrections of higher orders 
in M^HAM 2 ) were calculated. The singlet contributions 
appear for the first time at three-loop order. An example 
diagram is shown in Fig. [13] 

In Ref. (76) C\y has been determined at four loop or 
der in the effective Lagrangian approach with the help 
of the knowledge of the anomalous dimensions. This 
result was subsequently also obtained independently 
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Figure 13: Example three-loop singlet diagram. Solid lines denote top 
quarks, wavy lines are photons, twisted lines represent gluons and the 
dashed line is the Higgs boson. 

through the calculation of the four-loop order of the de¬ 
coupling function £ gy in Ref. f77j where also the corres¬ 
ponding five-loop contributions were determined, with 
the help of the anomalous dimensions. Some example 
diagrams at four-loop order are depicted in Fig.[f4| 



Figure 14: Example diagrams which illustrate different kind of dia¬ 
gram classes which have been determined in the heavy top-quark mass 
limit at four-loop order. Dotted lines represent massless quarks; all 
other lines are as defined in Fig.[l3] 


5. Conclusions 

The systematic investigation of four-loop tadpole in¬ 
tegrals started about a decade ago. In this article we 
have briefly touched the techniques which have been 
developed to perform the reduction to master integrals 
and to obtain results for the latter. Furthermore, we 
have described in some detail the most important ap¬ 
plications. Among them are four-loop corrections to 
the electroweak p parameter, the precise determination 


of charm and bottom quark masses and the decoupling 
constants in QCD. The latter have a close connection 
to Higgs boson production and decay into gluons and 
photons which has also been elaborated. 
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